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Abstract 

We develop a formalism for performing real space renormalization 
group transformations of the "decimation type" using perturbation 
theory. The type of transformations beyond d = 1 is nontrivial even 
for free theories. We check the formalism on solvable case of 0(N) 
symmetric Heisenberg chain. The transformation is particularly useful 
to study asymptotically free theories. Results for one class of such 
models, the d=2 O(N) symmetric a models (N > 3) for decimation 
with scale factor rj = 2 (when quarter of the points is left) are given 
as an example. 
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The renormalization group (RG) transformations is one of the most pow- 
erful and frequently used conceptual as well as practical tools in statistical 
physics and quantum field theory. While conceptually the idea of combining 
variables on neighbouring sites into complexes is very simple, in practise it 
almost always turns out to be rather complicated. 

Historically, usefulness of RG transformations was realized after the d=2 
Ising model on triangular lattice was very elegantly solved by Niemeijer and 
Van Leeuwen |J using block spinning. The nonlinear RG method they used 
was peculiar to that particular system and did not allow generalization to 
more complicated cases. For more general systems Wilson proposed to use 
the weak coupling perturbation theory in momentum space. This was first 
applied to scalar 4 models and subsequently to spin systems || and lattice 
gauge theories for thinning by factor 2 ||] . For complicated systems like these 
with local gauge symmetries some approximate methods were developed like 
Migdal-Kadanoff || approximation, variational RG ||, mean- field RG [|7| or 
block spinning using Monte Carlo numerical methods f| |9[ . However, unlike 
perturbation theory, these approximations are uncontrollable in a sense that 
it is not clear how to estimate errors. 

Exact RG transformations are generally not known (exceptions are dec- 
imations in spin chains d = 1 and mentioned above very special cases in 
d — 2). Moreover after one RG transformation the resulting action contains 
generically infinite number of interaction terms, and therefore one is forced 
to make an additional approximation dropping some of them (hopefully the 
less relevant ones). RG transformations are especially useful when applied 
repeatedly. This requires self similarity of the approximate effective action 
and is justified only around fixed points. Note however that accurate thin- 
ning of the lattice even just by factor r\ = 2 can greatly facilitate the study 
of a model by means of subsequent MC simulation. 

In this paper we perform the decimation for multidimensional free the- 
ories and then develop a systematic formalism to calculate perturbatively a 
decimated effective action for interacting models. [] 

There are several types of the RG transformations. The conceptually 
simplest one is the decimation or thinning of degrees of freedom in the con- 
figuration space. Some degrees of freedom located, for example, on sites with 

1 Method, based on the Schur formula, somewhat similar to ours, was developed for 
quasiperiodic systems on the octagonal lattice (Hi. 
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at least one odd coordinate are simply integrated out. 
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Fig. 1. Decimation: full circles belong to sublattice C*, 
while empty circles x G £> = C — C* denote the integrated out sites. 



Example is given on Fig. 1 on which spins at empty circle points are 
integrated out. 

Z = 2 e ~ A[ ^ x) '^ x)] = e ~ Adec[ ^ x)] (1) 

<p(X):XeC* cj>(x):xeC-C* 4>(X):X&C* 

Here and in what follows points of the coarse lattice are denoted by capital 
letters. The resulting effective action A dec contains generally interactions of 
any range. Here £ is the original d dimensional lattice (lattice spacing a), 
while C* is a sublattice with lattice spacing A = rja. In what follows A will 
set a scale: A — 1. Note that remaining variables are all the old variables. 
This is not the case for the so called block spin transformations. One defines 
a linear or a nonlinear combination of the variables on C, the block spin: 

<j ) (X) = f[X, ( / ) (x)} (2) 
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For example, for the O(N) classical spins S a one can define H] 



S a (X) = £ S a (x)/\ £ S a (x)\ (3) 

blockX blockX 

The choice of the combination is highly ambiguous and success of the trans- 
formation critically depends on it. The main problem is that it is extremely 
difficult in practise to perform such a transformation even perturbatively. 



The relations like eq.p) are very nonlinear and even singular 11 



Another type of RG transformations, used especially extensively in field 



theory, is the momentum space RG g, |T^] . One defines the momentum space 
variables 

(j){p) = l/(2n) d [ d d xe ipx (j)(x) (4) 



Now one performs integration over high frequence modes (strictly speak- 
ing chopping the Brillouin zone, but more often the approximate spherically 
symmetric momentum cutoff A is utilized [Q, |13|). This type of RG trans- 
formations, while convenient for the 4 model, turns out to be especially 
inconvenient for constrained systems like the O(N) symmetric a model. The 
reason is following. While generally in x - space the constraint are local, for 
example 

S a (x)S a (x) = 1 (5) 

in p - space it becomes a convolution. What does it mean now high frequency 
physical modes? The constraint mixes between low and high frequencies. 
Since most systems of interest belong to this class one has to circumvent the 
difficulty. One way is to solve the constraint and make the momentum space 
RG for physical quantities only. Then the mode integrated effective action 
contains generally "noncovariant terms". The original global symmetry is 
lost since the high frequency modes do not constitute an O(N) symmetric 
set. Problems are more acute with local gauge symmetries. In practice this 
type of thinning out of degrees of freedom is often used for the demonstration 
purposes only and very rarely the actual calculations. 

Decimation are extremely difficult to perform even in free theory in more 
then one dimension (see, for example |14|). It might sound surprising that 



something is difficult in free theory since all the integrals are gaussian and 
"doable in principle". Of course it is still a gaussian integral, but a very 
complicated one. Let us consider a free massless boson nearest neighbours 
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action 



A4>] 



Jd-2) 



J2<J)(x)d(x -y)(j){y) 



(6) 



where the lattice Laplacian is defined by = J2t=i [${ x ~ AO + $( x + lA ~ 25 (x)]. 

If one tries to integrate out a point 0(0, ...,0), the gaussian integral in- 
volves all its 2d nearest neighbours. 



J #(0)exp -^1 (2#(O) 2 -20(O)£(0(/i) + 0(-^) 



cxp 



2d 



V) 



(7) 



This is very simple. However when trying subsequently to integrate another 
point, say (2,0...,0), all the previous point's neighbours enter the gaussian 
integral and so on. The gaussian integration requires inverting increasingly 
larger matrices. Since we have to integrate out all the points not belonging to 
the sublattice, some other methods are required. An exception is the d = 1 
case. Here the size of the matrix does not grow: integration of a point leads 
just to interactions of the neighbouring remaining points. This is the reason 
why it is possible in many cases to explicitly find decimations in d — 1. 

Let us return with free massless boson on lattice eq.([j|). Effective action, 
after decimation with parameter r] generally has a form: 



A A [<j>{x)] = -i £ 4>(x)A(x - Y)4>(Y) 

xeR d 



(8) 



where bold letters denote sublattice functions. 

Using corresponding Fourier transforms we write propagators on the lat- 
tice and the sublattice: 



G{x) 



1 



(2vr) d J-n/a 



ir/a 



d d ke 



-ikxa 



4£ M sin 2 (a^/2) 



Z7T 



-iKX 



A(K) 



(9) 
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The two actions should give an equivalent correlators between two sub- 
lattice points: and X: G(r)X) = G(X). This leads to the following relation 
between the Fourier transforms: 



G(K) = 1/A(K) = j-*—; y e~ tKX d d ke lkX %, , , . (10) 

1 ^ ' K ' (2it) d j? J-*f* E M 4sin 2 (a^/2) 1 ' 

Summation over X results in sum over 6 functions 

G <*> - C* k JL m ~ K + 2m) ^^w) 

which are used to perform the integrations over k: 

v a 2 
G(K) = Y ^— — - (12) 

Limits of summation follow from the different sizes of Brillouin zone for two 
lattices. In d = 1 we recover the original form of the action, while for d > 1, 
the sum over one of the variables, n\ can be explicitly done |L5[], but the 
remaining summations should be done numerically. In particular for d = 2 
the propagator is: 

= I f sinh(^)csch(q) 

V £i 2 (-1 + cosh^a) + 2sin 2 (f ) X 



,2 



where 



2 K 2 + 2im 2 . 



a = arccosh(l + 2sin ( )). 

2r] 



For large rj the euclidean invariance is restored, G(Ki, K 2 ) = G(y K\ + K 2 ) 
and numerical calculations show that it can be fitted by 

G{K) = ^- + — log(?7) + 0.04876 + 0.003022 AT 2 + O [{K 2 ) 2 ] (14) 
K 1 2tt 

with an accuracy better 1 percent in all the Brillouin zone. Note that the 
decimated propagator even for large r) does not coincides with the naive con- 
tinuum limit 1/K 2 . The contact constant term with logarithmic dependence 
of rj is typical for d = 2 and is nothing else but the bubble integral. The 
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polynomial coefficients are very small and almost coincide with the Loran ex- 
pansion of the propagator around K — 0. For finite 77 the symmetry remains 
of course just the discrete subgroup of the rotations. In higher dimensions 
similar expressions can be written. The procedure can be extended also to 
free fermion fields. 

Now we briefly outline a systematic procedure for decimation-type RG 
transformations in interacting models. For concreteness we discuss the lattice 
4 model 



;(V0,) 2 + 



m 2 a 2 



+ 



Xa 2 



4! 



where V(fi x = {(fi x - (fi x -i). 

In the process of decimation, the fields (fix on the sublattice (X 
will be treated as "external sources". 



(15) 



G £*) 



O- 
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Fig. 2. Real space RG (decimation) propagators (a, b, c) and vertices (d, e) 
for 4 model. Full circles belong to sublattice (external fields), 
while empty denote "internal fields" (x G £ — £*). 



Action can be divided into three parts: "classical action" of the external 
sources (fix (Fig. 2a, 2d), 

-—(fix + a d 2 d (fi 2 x + 
x \ A 



A d [<fix] = E 



4! 
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= £ <PxA XY <p Y + AiM, (i6) 

Z X,Y 

cross-term (Fig. 2b) 

A-y[<f> x , <f>x] = -a {d - 2) E 0xA(X, x)<f> x = -a^ £ <j> x B xX <t> x (17) 

x,X 

with "external legs" 

A(X, x) = + (18) 

and an internal part for which all the vertices belong to "D = C — C* (" deco- 
rated" model) (Fig. 2c, 2e): 

AW = E ^ + ir^^ (19) 

" x,y ^ - x 

Note that, unlike momentum space RG, here external fields <px are coupled to 
internal part only via derivative couplings like off-diagonal part of propagator 
(Fig. 2b) or derivative interaction in nonlinear cr-model. All the local vertices 
will completely decouple into internal (Fig.2e) and external (Fig. 2d). 

Integration out all the fields 4> x will lead to an effective action for the 
fields 4> x on sublattice of the form: 

A ^x]= E ^# (2n HAV.,A^0 Xl ...0 X2n , (20) 

where the coefficient functions H < - 2n \X 1 , ..,X 2n ) are sums of all the connected 
contributions with 2n ends. These connected functions does not degenerate 
into one particle irreducible. 

To integrate over field <p x perturbatively, we need to find its propagator 
which is matrix inverse to D. At the same time, all we explicitly have is 
the full Laplacian A and original propagator G = A -1 . Moreover, since 
V does not constitute a sublattice, it is impossible to make use of Fourier 
analysis on it. Therefore it is useful to represent all the summations over D 
via summations over C and C*. This can be done using following algebraic 
trick. Matrix A as well as matrix G can be decomposed into following blocks 

A=(£ I). (2D 
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°-(? 5). ^ 

Infinite-dimensional matrices A, B, D defined by the quadratic part of action 



eqs. (jig , |18| , |19|) (Fig. 2a, 2b, 2c). The matrices b,d are the usual propa- 
gator matrices (Fig. 2b, 2c) and a is the propagator between the points of 
sublattice, that is an expression inverse to G(X). Now we can invert D using 
the fact that matrices A and G are inverse to each other. This implies a set 
of algebraic relations for their submatrices: 

Aa + Bb l = 1, Ab + Bd = 0, B l a + Db l = 0, B% + Dd = 1 

and, after straightforward transformations we obtain an expression for the 
"internal" propagator: 

IT 1 =d- b l ab (23) 
or, returning to previous notations, 

D xy - Gxy - GxxGxyGyy (24) 
X,Y 



Here and in what follows bold line denotes free inverse decimated propa- 
gator Gxy, while thin line corresponds to propagator of the original theory 
G xy . Using this representation of D~y as an internal line, we can now begin 
to build the perturbation theory. We would like to stress that using of D~y 
eq. (|24l) enables us to extend this summation over all the original lattice. In- 
deed one can see that this expression is equal to zero when at least one of 
the points x, y in belongs to C* 

Calculation of an n-point function H^ n \X 1 , .., X n ) in effective action 
eq. (PTJ|) for the coarse-grained field <px will be as follows. 

All connected diagrams with n end points Xi,..,X n , V vertices and / 
internal lines in real space are drawn with following components: 

a) all vertices are situated at the points x E £ and to every vertex at 
point Xi corresponds summation ^2 Xi ] 

b) Axx are assigned to external ends X and 

c) the internal lines D~y are represented via eq. (|24|) . This representation 
splits each diagram into 2 1 subdiagrams and each of these subdiagrams should 
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be calculated separately. This calculation includes summation over all the 
internal points Xi on the fine grained lattice and over the internal sublattice 
points Y [end points of the inverse decimated propagator G" 1 , see eq.(P^)]. 
Finally the "classical contribution" to coefficient function should be added. 

For the sake of simplicity let us discuss calculations of the decimation 
diagrams on the concrete example. Namely, we will consider one of the 
contributions to four point function in 4 model (Fig. 3). 

With using of the propagator D~y, original diagram splits into five sub- 
diagrams (Fig. 3a, 3b, 3c, 3d, 3e). Typical subdiagram here (for instance, 
sub diagram d) can be written as 

Hf\Xi, ..,X 4 ) = 

^2 Axi yiG yi y 1 G y Iy 2 Gy 2 zAx 2 y 2 G y2 y a GyIy a Gy4, z^X 3 y 3 G X3 2^X4 yiG Xi z . 

(25) 

As usual, in practical calculations it is very convenient to employ the Fourier 
transformed functions at intermediate steps. In this way we will operate with 
vertices, legs 

A(fc) =2a d - 2 ]Tcos(£; M a), (26) 

propagators G(k) and inverse decimated propagator G~ 1 (ii') eq. fljjD . How- 
ever, it is not convenient to perform Fourier transform of expression eq.(|2"5|) 
immediately, because it contains functions defined on different lattices. 
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(a) 




Instead, we can use the fact that this expression breaks into blocks where 
all internal points lie on the fine grained lattice C and sublattice points enter 
only as ends. These blocks are connected by G~ 1 (Y i: Yj). In terms of such 
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an " /^-connected" parts the diagram eq.(p5|) has the form: 

H^(X U ..,X 4 ) = £ F 1 (X 1 , Y 1 )F 2 (X 2 , Y 2 )F 3 (X 3 , X 4 , Y 3 , Y A )G^ Ya G^ Y 

Yt.,..Y 4 

(27) 

with blocks 

F 1 (X 1 , Y 1 )=J2^x iyi G yiYl , (28) 

yi 

F 2 (X 2 , Y 2 ) = J2^x 2y2 G y2 Y 3 (29) 

y-2 

and 

-^3(^3, -^4) ^3> Y4) = Gy 2 zGy 4 zAx 3 y3 G X3 z Ax 4 y 4 G Xi z . (30) 

The most economic way to proceed is to Fourier transform each block in 
eq. (|27|) with respect to the fine grained lattice variables only. This will result 
in diagram as function on sublattice, and then we can apply to it the 
Fourier transformation rules for the sublattice. First stage does not differ 
from the usual perturbative lattice calculations. Resulting functions are, for 
example, as 

F X {X 1 , n) = r ° e i(Xl - yi)fc A(A;)G(A;). (31) 

J —it /a 

This, however is not the case when the remaining Fourier transforms and 
summations over internal sublattice points are performed. Due to the dif- 
ference between the Brillouin zones result will be sums rather then simply 
products. For instance, for the block F\ we will obtain: 

all n M =l 

Such an expressions we will call "decimated" and denote by [|"expr"|]: 

[\f(K,L,. ..){}= £ f(K + 2nn K ,L + 27in L ,...) (33) 

(n K ^,n Lti ,..)=l 

Notice that the decimated function does not possess original translation in- 
variance with the period 27r/a, but instead period is 27i/(r]a). This corre- 
sponds to the nature of decimation as transformation £—>■£*. Eventually, 



E 



,(d-2) 



COS 



+ 2im p 



_ 77 2 £„.4sin 



2/ Ku+2-nnn 



(32) 
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diagram takes the following form: 

Hf\K,P,Q) = [|F 1 (if)|]G- 1 (if)[|F 1 (P)|]G- 1 (P)[|Fs(lf I P,g,-iir-P-g)|]. 

(34) 

The described algorithm is applicable not only to the tree level contri- 
butions. When the representation of the propagator eq. (|24]) is used, some 
of the loop diagrams will contain G -1 . If this is the case, one should calcu- 
late the corresponding decimated functions, and only then perform the loop 
integration. 

We applied the formalism of the real space RG to O(N) symmetric non- 
linear sigma model. In d dimensions, this model is described by action 

A = -—^2s a ns a , ( 35 ) 

where S a (a = 1, .., N) is O(N) vector normalized on unity, S 2 = 1, g is the 
coupling constant. The most general covariant effective action with up to 
four lattice derivatives is: 



1 -S a us a + c 5 (as a ) 2 + {— + c 6 ) £(<V^) 2 



+c 7 (S a as a f + c 8 SC^flf S a ) 2 + c 9 W Jt^-JL S a .^-^S a 

H nv \ 1 1 ) 

(36) 

Note that Symanzik- improved action |l6j contains O(N) noncovariant terms. 
These can be covariantized |17j introducing source terms like JOS, J 2 . In 
real space decimation these obviously do not appear. 

Decimation technique is applied to the unconstrained variables, while 
both original and effective actions eqs. (|35|j36| ) are written in terms of con- 
strained, covariant fields S a . To build a perturbation theory, one can solve 
the constraint: 

S a = (v^^V, y/l-gaWfr*)*), i = l,..,N-l, 

and expand these actions in terms of "pions" 7r\ Then to the fourth order 
in 7T l and up to fourth derivatives, the effective action is: 
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with the quadratic part 



E 

X 



Z4 A« 



(37) 



and the quartic part 

^ (4) = -sE 



+ c 6 ) E(^ 2 ) 2 



24 



+c 7 (ttW) 2 + c 8 E(^^ 5 >T + c 9 E 



d+ . d +d + 



(38) 



To reconstruct the coefficients C5,..,Cg, we perform the perturbative decima- 
tion of the original theory, obtaining the decimated coefficient functions H^ n > . 
Then the coefficients can be identified by simple comparison. We applied the 
method to d = 2 a model for rj = 2. Comparing the quadratic parts at tree 
level, we find c 5 , c 6 : 

c 5 = 5/32; c 6 = 1/96. (39) 

There are three unknown coefficients in the quartic part of the effective action 
eq.(|38p. To determine them, we calculated the quartic decimated function 
if ( 4 Mor three different momenta configurations, thus obtaining a system of 
the linear equations for C7, Cs, C9 with the solution: 



c 7 = c 8 = 0, c 9 



-5/64. 



(40) 



The vanishing of C7 and c§ is rather surprising. We do not see any obvious 
reason for this. 

To see how formalism works we applied it to the solvable d = 1 O(N) 
symmetric a model at both tree and one loop level. In this case there is only 
one quartic term, |f X)x( 7r x'-' 7r x) 2 m the effective action. We calculated two 
and four point functions. These perturbative RG results were compared to 
expansion of exact solution for 77 = 2 decimated effective action of Heisenberg 
chain [18]: 

-A e " = (l/(2^)-(iV-2)/4)(^^ +1 -l)-(l/(16^)-(iV-2)/16)(^^ +1 -l) 2 +. 

(41) 

and are in complete agreement. 
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We do not present any technique details here. Only note, that at the 
first sight the decimated coefficient functions should behave as 1/P 2 
and as 1/P 4 . In fact, however, due to the decimation procedure all 
the negative powers of momenta and the constant terms cancel. This exact 
cancellations provide an additional consistency check. 

To summarize we found a systematic way to perform RG of the decimation 
type in d > 1 perturbatively. Here we would like to discuss some of its uses. 
The formalism we propose here preserves all the local relations including 
constraints. This in turn means that effective (coarse grained) theory will 
obey exactly the same local constraints as did original, and that no non- 
covariant terms will appear in the effective action. 

Another, compared to others 0, |5], [|, useful feature of proposed for- 
malism is its perturbative character. This can provide us with systematic 
method of calculations in asymptotically free models and, what is even more 
essential, with a way to do controllable approximations. Hopefully, this side 
of proposed formalism will make it applicable in situations when such a con- 
trol is essential, as in the recently proposed double strong-weak expansion 
approach |L8||. It turns out that in asymptotically free theories there exist a 
region in the parameter space in which both strong and weak coupling expan- 
sions are valid at the same time. Namely both the practical weak coupling 
a(g) = const g and strong coupling (3(g) = const jg expansion parame- 
ters there are reasonably small. The "loop factors" l/(4ir) 2 in the practical 
weak coupling expansion parameter a(g) are partly responsible for this. In 
this scheme high frequency modes are integrated out perturbatively and the 
resulting effective action treated using strong coupling expansion. The sym- 
metry preserving perturbative decimation technique is the most suitable tool 
for the first part of such calculations. 

We would like to stress also that the method described here unlike most 
of the other decimation (and exact RG in general) techniques, enables us 
to perform decimations not restricted to the simplest case rj — 2 only. The 
i] > 2 calculation just takes a bit more computer time. 

Discussions with B. Hu and H. Miller are greatly appreciated. This work 
was supported by National Science Council of ROC, grants NSC-83-0208-M- 
001-011 (B.R.) and NSC-83-0208-M-001-015 (V.K.). 
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